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Abstract 

We present a partonic Monte Carlo event generator that can be used for calculating the group 
independent kinematical functions of any infrared safe four-jet observable at next-to-leading 
order accuracy. As an example, we calculate the differential distribution of the IIi and 114 Fox- 
Wolfram moments. We find large K factors (K > 2). The effect of the radiative correction is 
to increase the overall normalization, but not to reduce the renormalization scale dependence 
significantly. 
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We present a partonic Monte Carlo event generator that can be used for calculating the group independent 
kinematical functions of any infrared safe four-jet observable at next-to-leading order accuracy. As an example, we 
calculate the differential distribution of the fli and II4 Fox- Wolfram moments. We find large K factors (K > 2). 
The effect of the radiative correction is to increase the overall normalization, but not to reduce the renormalization 
scale dependence significantly. 



1. INTRODUCTION 

The analysis of four-jet production in e+e~ an- 
nihilation is an important tool to learn about the 
the basic properties of Quantum Chromodynam- 
ics (QCD), the theory of strong interactions. So 
far the four-jet data were used mainly for color 
factor measurements because the theoreti- 
cal prediction of perturbative QCD has not been 
known at next-to-leading order accuracy that is 
needed for precision as measurements. In this 
contribution we report the construction of a par- 
tonic event generator that can be used for calcu- 
lating the radiative corrections to the group inde- 
pendent kinematical functions of any infrared safe 
four-jet observables in electron positron annihila- 
tion. These results make possible the simultane- 
ous precision measurement of the strong coupling 
and the color charge factors using LEP or SLC 
data. As an example, we present the differential 
distributions of two event shape variables that are 
non-trivial for four-jet like events — the Hi and 
H4 Fox- Wolfram moments 
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2. THE METHOD 

The higher order correction to the leading order 
partonic cross section is a sum of two terms, the 
real and virtual corrections: 

a^LO ^ J d^NLO ^ ^ ^^R + da^ , (1) 

where in case of four-jet observables, the integrals 
are over the five- and four-particle phase space 
respectively. These two integrals are both diver- 
gent in d = 4 space-time dimensions, however, 
their sum is finite for infrared safe physical quan- 
tities. In recent years several general methods 
have been developed for exposing the cancella- 
tion of infrared divergences directly at the inte- 
grand level ||-|]. We use a modified version of the 
dipole formalism of Catani and Seymour jsj.The 
formal result of this cancellation is that the next- 
to-leading order correction is a sum of two finite 
integrals, 

aNLO^^d^NLO_^^d^NLO^ (2) 

where the first term is an integral over the avail- 
able five-parton phase space (as defined by the jet 
observable) and the second one is an integral over 
the four-parton phase space. The distinct feature 
of this formalism as compared to other cancella- 
tion methods is that a single subtraction term is 
used for the regularization of the real cross sec- 
tion. This subtraction term provides a smooth 
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approximation of the real cross section in all of 
its singular limits (soft and coUinear regions), re- 
sulting in a well-converging partonic Monte Carlo 
program. 

The main ingredients of the calculation are the 
four-parton next-to-leading order and five-parton 
Born level squared matrix elements. The helic- 
ity amplitudes from which the latter can be con- 
structed have been known for almost a decade 
Recently, there was important development 
in the calculation of the virtual corrections for the 
processes e^e~ — > qqQQ and qqgg. On one hand 
Campbell, Glover and Miller made FORTRAN 
programs that calculates the next-to-leading or- 
der corrections to the four parton processes via 
the production of an s-channel virtual photon 
publicly available |Q, while on the other, the 
new techniques developed by Bern, Dixon and 
Kosower in the calculation of one-loop multipar- 
ton amplitudes Q made possible the derivation of 
explicit analytic expressions for the helicity am- 
plitudes of the e+e~ Z", 7* qqQQ and qqgg 
processes These results — and so ours — 

are valid in the limit when all quark and lepton 
massess are set to zero. We use the helicity am- 
plitudes of refs. [9IJ10I] for the loop corrections. 

3. GENERAL STRUCTURE OF THE 
CROSS SECTION 

Once the integrations in eq. are carried out, 
the next-to-leading order differential cross section 
for a four-jet observable O4 takes the form 
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In eq. (^) (Tq denotes the Born cross section for 



the process e+e 
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with the normalization = 1/2 in Tr(T"Tt'') = 
TuS""^, s is the total cm. energy squared, n is 
the renormalization scale, while Bq^ and C04 are 
scale independent functions, i3o4 is the Born ap- 
proximation and is the radiative correction. 
The Born approximation and the higher order 



correction are linear and quadratic forms of ra- 
tios of eigenvalues of the Casimir operators of the 
underlying gauge group 



B04 = Bq + Bj;X + By-y , 
and 

C04 ^Co+ C^x + Cyy + C,z 
+ Cxx X -\- Cxy xy ~\- Cyy y 



(4) 
(5) 



The X and y parameters are ratios of the quad- 
ratic Casimirs, x — Ca/Cf and y = Tn/Cp, 
while z is related to the square of a cubic Casimir, 
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C3 = ^ Tr(T°T''r^'=)Tr(r^^T''r") 



(6) 



,,b,c=l 

- g3 



4. RESULTS 

Dixon and Signer obtained the first ever com- 
plete results for four-jet observables at next-to- 
leading order accuracy [l^ , which were four-jet 
rates for three different clustering algorithms: the 
Durham Hi, the Geneva and the EO ^ 
schemes. In order to compare the two programs, 
we have also calculated these quantities at the 
same values of the parameters and found very 
good agreement 

As a new result, we present the next-to- leading 
order prediction for Fox- Wolfram moments Hi 
and 114 (H . These observables were constructed in 
such a way that they vanish for coplanar events, 
thus have non-trivial value for four-jet like events. 
The variable Hi is defined as 
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1 = 
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where pi is the three-momentum of particle i, is 
the unit vector along the momentum pi, and the 
sum runs over all final state particles in an event. 
The variable 114 has more complicated definition: 
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Figure 1. Next-to-leading order corrections to the 
group independent kinematical functions of the 
event shape variable Hi. 



Figure 2. Next- to-leading order QCD prediction 
of the event shape variable IIi. 



We present histograms for the various group 

independent kinematical correction functions Ci 
for Hi in Figure 1. We find that the correction 
functions Cq and are large and positive, while 
the functions Cy and Cyy are large and negative. 
We have not shown the contribution, that is 
negligible. 

Wc take SU(3) as underlying gauge group, and 
obtain the next-to-leading order prediction for 
the moments using the QCD values x = 9/4 
and y = 3/8. The results — obtained for five 
light quark flavors at the Z° peak with Mz = 
91.187GeV, Tz = 2.49GeV, sin^ Ow = 0.23 and 
as{Mz) = 0.118 — are plotted in Figure 2 and 
3. The light grey bands indicate the renormal- 
ization scale dependence in the range 0.1 < a;^ = 
iil\fs < 2 at leading order, while the dark bands 
show that at next-to-leading order. We see that 
the effect of the higher order correction is to in- 
crease the overall normalization, but no signifi- 
cant scale-dependence reduction occurs. 

In order to see the renormalization scale depen- 
dence in more details, we define the average value 
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Figure 3. Next-to- leading order QCD prediction 
of the event shape variable 114. 
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Figure 4. Renormalization scale dependence of ' 
the average value of the moment Hi. 

Figure 5. K factors of the Hi and 114 moments. 



of these shape variables as 

<04>=- f' dO^O^^ , (9) 

and study the dependence of the average value of 
the Hi variable on the scale in Figure 4. We see 
that there remains substantial scale dependence 
at next-to-leading order showing that the uncal- 
culatcd higher order corrections are presumably 
large. The feature is similar for the 114 variable, 
but the residual scale dependence is even larger. 

The same conchision is drawn if we look at the 
dependence of the K factors on the observables 
as depicted in Figure 5. In case of the IIi pa- 
rameter the K factor is slightly above two for the 
whole range, while for 114 it is even larger. This 
suggests that 114 cannot be reliable calculated in 
perturbation theory. 

5. CONCLUSION 

In this talk we presented for the first time a 
next-to-leading order calculation of the differen- 
tial cross section of two four jet shape variables, 
the Hi and 114 Fox- Wolfram moments. We gave 
explicit results both for the full next-to-leading 
order cross sections and for the group indepen- 
dent kinematical functions of the radiative correc- 
tions. The effect of the radiative correction was 
to increase the overall normalization, but not to 
reduce the renormalization scale dependence sig- 



nificantly. 

These results were produced by a partonic 
Monte Carlo program DEBRECEN that can be 
used for the calculation of QCD radiative correc- 
tions to the differential cross section of any kind 
of infrared safe four-jet observable in electron- 
positron annihilation. 

We thank L. Dixon for providing us the Math- 
ematica files of the one-loop four-parton helicity 
amplitudes. 
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6. DISCUSSIONS 

A. P. Contogouris, Universtity of Athens 
To produce your K factors you need a jet algo- 
rithm; however, there are more than one algo- 
rithms. Do I understand correctly that your K 
factors much depend on it? 

Z. Trocsanyi 
The K factors do depend on the observable. This 
is also the case for three-jet observables. The 
smaller K factor the more reliable the quantity 
can be calculated in perturbation theory. From 
this point of view the Durham jet clustering al- 
gorithm (K'^1.6) is a better observable than the 
event shape variables. 

A. P. Contogouris 
Suppose you compare some of your K factors with 
those for three-jet production in e^e^ annihila- 
tion, corresponding to equivalent observables (e.g. 
shape variable). How do your K factors compare? 
Also, how do they compare with Drell- Yan lep- 
ton pari production, one of the first processes for 
which K- factors were calculated? 

Z. Trocsanyi 
The K factor for the C parameter for instance, 
is about 1.5, i.e. much smaller than in our case. 
The Drell- Yan K factor is 1.8-2, which is again 
quite large, but somwhat smaller than our values. 



